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ABSTRACT 

In the context of a relativistic hot spot model, we investigate different physical 
mechanisms to explain the behavior of quasi-periodic oscillations (QPOs) from 
accreting black holes. The locations and amplitudes of the QPO peaks are de- 
termined by the ray-tracing calculations presented in Schnittman & Bertschinger 
(2004a): the black hole mass and angular momentum give the geodesic coor- 
dinate frequencies, while the disk inclination and the hot spot size, shape, and 
overbrightness give the amplitudes of the different peaks. In this paper addi- 
tional features are added to the existing model to explain the broadening of the 
QPO peaks as well as the damping of higher frequency harmonics in the power 
spectrum. We present a number of analytic results that closely agree with more 
detailed numerical calculations. Four primary pieces are developed: the addition 
of multiple hot spots with random phases, a finite width in the distribution of 
geodesic orbits, Poisson sampling of the detected photons, and the scattering of 
photons from the hot spot through a corona of hot electrons around the black 
hole. Finally, the complete model is used to fit the observed power spectra of 
both type A and type B QPOs seen in XTE J1550-564, giving confidence limits 
on each of the model parameters. 

Subject headings: black hole physics - accretion disks - X-rays:binaries 



1. INTRODUCTION 

One of the most exciting results from the Rossi X-Ray Timing Explorer (RXTE) was 
the discovery of high frequency quasi-periodic oscillations (HFQPOs) from neutron star and 
black hole binaries (Strohmayer et al. 1996; Strohmayer 2001; Lamb 2002). For black hole 
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systems, these HFQPOs are observed repeatedly at more or less constant frequencies, and 
in a few cases with integer ratios (Remillard et al. 2002; Homan et al. 2004; Remillard et al. 
2004). These discoveries give the exciting prospect of determining a black hole's mass and 
spin, as well as testing general relativity in the strong-field regime. 

To understand these observations more quantitatively, we have developed a ray-tracing 
code to model the X-ray light curve from a collection of "hot spots," small regions of excess 
emission moving on geodesic orbits (Schnittman & Bertschinger 2004a,b). This hot spot 
model is motivated by the similarity between the QPO frequencies and the black hole co- 
ordinate frequencies near the inner-most stable circular orbit (ISCO) (Stella & Vietri 1998, 
1999), as well as the suggestion of a resonance leading to integer commensurabilities be- 
tween these coordinate frequencies (Abramowicz & Kluzniak 2001, 2003). Stella & Vietri 
(1999) investigated primarily the QPO frequency pairs found in LMXBs with a neutron star 
accretor, but their approach can be applied to black hole systems as well. 

The basic geodesic hot spot model is characterized by the black hole mass and spin, the 
disk inclination angle, and the hot spot size, shape, and overbrightness. Motivated by the 
3:2 frequency commensurabilites observed in QPOs from XTE J1550-564, GRO J1655-40, 
and H1743-322 (Remillard et al. 2002; Homan et al. 2004; Remillard et al. 2004), we pick a 
radius for the geodesic orbits such that the coordinate frequencies and v r will have a 3:1 
ratio. Thus the 3:2 commensurability is interpreted as the fundamental orbital frequency 
and its beat mode with the radial frequency at — v r . Together with the assumption that 
the concurrent low frequency QPOs are due to Lense-Thirring precession at this same radius 
(Stella et al. 1999), the location of the QPO peaks uniquely determines the black hole mass 
and spin. Relaxing the low frequency QPO criterion leaves a one-dimensional degeneracy 
in the mass-spin parameter space which can be broken by an independent determination of 
the binary system's inclination and radial velocity measurements of the low-mass companion 
star [e. g., Orosz et al. (2002, 2004)]. 

Given the black hole mass, spin, inclination, and the radius of the geodesic orbit, the 
parameters of the hot spot model (i. e. the hot spot size, shape, and overbrightness, and the 
orbital eccentricity) are determined by fitting to the amplitudes of the peaks in the observed 
power spectrum. However, the model as described so far produces a perfectly periodic X-ray 
light curve as the single hot spot orbits the black hole indefinitely. Such a periodic light 
curve will give a power spectrum composed solely of delta-function peaks (Schnittman & 
Bertschinger 2004a), unlike the broad features in the observations. 

In this paper we introduce two simple physical models to account for this broadening 
of the QPO peaks. The models are based on analytic results, then tested and confirmed 
by comparison with the three-dimensional ray-tracing calculations of multiple hot spots 
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(Schnittman & Bertschinger 2004b). These results draw repeatedly from the convolution 
theorem of complex Fourier analysis, an essential tool for relating the behavior of the physical 
system in the time domain with the more easily measured power spectrum in the frequency 
domain. We find the power spectrum can be accurately modeled by a superposition of 
Lorentzian peaks, consistent with the standard analysis of QPO data from neutron stars and 
black holes (Nowak 2000; Belloni et al. 2002). Many of the methods and results presented 
here are equally valid for other QPO models such as diskoseismology (Wagoner, Silbergleit, & 
Ortega- Rodriguez 2001), vertically-integrated disk oscillations (Rezzolla et al. 2003), toroidal 
perturbations (Lee & Ramirez- Ruiz 2002; Lee et al. 2004), and magnetic resonances (Wang 
et al. 2003). 

In Section 2, we derive the effect of summing the light curves from multiple hot spots 
with random phases and different lifetimes to give a Lorentz-broadened peak in the power 
spectrum. Section 3 shows how a finite width in the radii of the geodesic orbits produces a 
corresponding broadening of the QPO peaks. In Section 4 we develop two simple models for 
detector statistics and photon scattering, both of which affect other features of the power 
spectrum such as the continuum noise and the damping of high frequency harmonics, but 
do not contribute to the broadening of the QPO peaks. Finally, all the pieces of the model 
are brought together in Section 5 and used to interpret the power spectra from a number 
of observations of XTE J1550-564. We conclude with a discussion of open questions for the 
hot spot model as well as directions for future work. 



2. PEAK BROADENING FROM RANDOM PHASES 

In this section we derive the shape of a QPO peak in Fourier space broadened by the 
summation of multiple periodic functions combined with random phases. There are many 
different accepted conventions for discrete and continuous Fourier transforms (Press et al. 
1997), so we begin by defining the forward- and reverse-transforms between the time and 
frequency domains (t and v). For a Fourier pair f(t) and F(u), 

Fj = i, X; /^"W - F(u) = -"r P f(t)e-™*dt (1) 

and 

f k =J2 F jf FW N - -> f(t) = T S / F{y)<***d», (2) 

where f(t) is defined on the time interval [0,7/] and v N = 1/(2 At) is the Nyquist frequency 
for a sampling rate At = Tf/N 8 . With this convention, f(t) and F(v) conveniently have the 
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same units and Parseval's theorem takes the form 

Tf rv N 

f 2 (t)dt = T 2 F\v)dv. (3) 

J —vn 

For such a time series f(t), the power spectrum is defined as F 2 (u), the squared amplitude 
of the Fourier transform. 

Consider a purely sinusoidal function 

f(t) = A sm(2nu t + (/)), (4) 

where is some constant phase. If there are an integer number of complete oscillations 
within the time Tf or in the limit of Tf — > oo, the Fourier transform of f(t) will be 

£ e i(*-*/2) v = j/q 
F(u) = { A e ~i{4>-*/2) v = - Vq . (5) 
otherwise 

If we then truncate the function fit) by multiplying it with a boxcar window function w(t) 
of length AT, the convolution theorem gives the transform of the resulting function g(t): 

g{t) = f(t)w(t) & G(u) = (F * W){v). (6) 

In the case where the window function is longer than a single period and short compared to 
the total sampling time (1/Vo < AT <^ Tf), the convolved power G 2 (v) can be approximated 
by 

2 ^ A 2 sin 2 [7r(z/±z/ )AT] 
G [V) ~ 47? tt 2 (u±u ) 2 ■ (7) 

For f(t) real, F(—u) = F*(y) and since we are primarily concerned with the power 
spectrum F 2 (—u) = F 2 (u), we will consider only positive frequencies in the analysis below 
(unless explicitly stated otherwise). Of course, when calculating the actual observable power 
in a signal, both positive and negative frequencies must be included. 

All the information about the phase of f(t) and the location in time of the window 
function is contained in the complex phase of the function G(u). This phase information is 
important when considering the total power contributed by a collection of signals, each with 
a different time window and random phase. When summing a series of complex functions 
with random phase, the total amplitude adds in quadrature as in a two-dimensional random 
walk. Therefore combining iV different segments of f(t), each of length AT and random 0, 
gives a Fourier transform with amplitude \^N\G(u)\, and thus the net power spectrum is 
NG 2 (v). 
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The result in equation (7) is valid only if every segment of f(t) has the exact same 
sampling length AT and frequency u . This is not necessarily the case in the context of 
the hot spot model, in which small regions of excess emission move along geodesic orbits, 
modulating the flux periodically as photons are relativistically beamed toward the observer. 
Three-dimensional MHD simulations of accretion disks suggest that in general such hot spots 
are continually being formed and destroyed with random phases, with a range of lifetimes 
and orbital frequencies (Hawley & Krolik 2001; De Villiers, Hawley, & Krolik 2003). 

For now, we consider the contribution from identical hot spots, assuming that each one 
forms around the same radius with similar size and overbrightness and survives for some 
finite time before being destroyed. Over this lifetime, the hot spot produces a coherent 
periodic light curve sampled by a window function of duration T. Analogous to radioactive 
decay processes, we assume that during each time step dt, the probability of the hot spot 
dissolving is dt/Ti, where TJ is the characteristic lifetime of the hot spots. In this case, the 
differential probability distribution of lifetimes T for coherent segments is 

P(T)dT= d ^e-^. (8) 

Over a sample time Tf ^> 7], if there are an average of iV S p t hot spots in existence at any 
time, the number of hot spots formed with a lifetime between T and T + dT is given by 

dN(T) = N spot ^e~ T ^dT. (9) 

Assuming for the time being that each coherent section of the light curve is given by 
the sinusoidal f(t) used above, we can sum all the individual segments to give the total light 
curve I(t) with corresponding power spectrum 

f oo 



/» = / G 2 (u,T)dN(T) 
Jo 



- ^\2nT f ) Jo ^V^o) 2 Tf 
T 1 

= 2iV t A 2 ^ 1 + 47r27 , 2(i/ _^ )2 . (10) 

Hence we find the shape of the resulting spectrum is a Lorentzian peaked around vq with 
characteristic width 

Ai/= — — . (11) 



If this model is a qualitatively accurate description of how hot spots form and dissolve 
in the disk, one immediate conclusion is that the oscillator quality factor Q can be fairly 
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high even for relatively short coherence times: 

Q ^fwm = ^- (12 > 

For example, if every hot spot has a lifetime of exactly four orbits (TJ = 4/z/ ), the central 
peak of the power spectrum G 2 (u, 7}) has coherence Q = 4.5, about what one would expect 
from a first-order estimate. However, after integrating over the exponential lifetime distri- 
bution to get the Lorentzian profile of equation (10), the resulting quality factor is Q = 12.6, 
roughly a factor of three higher. Remillard et al. (2002) observe quality values of Q ~ 5 — 10 
for the HFQPOs seen in XTE J1550-564, corresponding to typical hot spot lifetimes of only 
2-3 orbits. 

In addition to the boxcar window, we have also tried other sampling functions to model 
the lifetime and evolution of the hot spots. Since the boxcar window represents an instan- 
taneous formation and subsequent destruction mechanism, the resulting power spectrum 
contains significant power at high frequencies, a general property of discontinuous functions. 

A smoother, Gaussian window function in time gives a Gaussian profile in frequency 
space: 

W W = ex P (^) = exp (2^) (13) 

where again the characteristic width is given by Av = 1/ (27rT) . After integrating over the 
same distribution of lifetimes dN(T) as above, we get the power spectrum 



j» = 4vriV spot ,4 2 -^ 3 



V^r(l + 2z 2 )erfc(z)e z2 - 2z , (14) 



where we have defined 

z = ^r, 7- (15) 

47rT,(i/ - i/ ) V ; 

For large z (near the peak at v — u Q ), equation (14) can be approximated by the narrow 

Lorentzian 

P(u) « 4nN svot A 2 — -L 777. (16) 

V ; spot T f 1 + 48tt 2 T z 2 (z/- z/ ) 2 

As with the boxcar window, the exponential lifetime distribution has the effect of nar- 
rowing the peak of the net power spectrum compared with that of a single segment of the 
light curve with length 7}. This smaller width can be understood by considering the dis- 
tribution of hot spot lifetimes [eqn. (9)] and their relative contribution to the total power 
spectrum [eqns. (7, 13)]. While there are actually more segments with individual lifetimes 
shorter than TJ, the few long-lived segments of the light curves add significantly more weight 
to the QPO peaks since W{y = 0,T) oc T while Av oc 1/T. 
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In fact, for any set of self-similar sampling functions w(t, T) — w(t/T), the corresponding 
power spectra W 2 (u,T) can be approximated near v = as a Lorentzian: 

w >' T)x TfTTWv' (17) 

with p a dimensionless constant over the set of w(t, T). The characteristic width of W 2 (u, T) 
is thus defined as l/(/3T). Integrating over the lifetime distribution dN(T), the net power 
function is given by 

/>W>,) 1 + 12/ 4 ( ,_, o)2 . (18) 

We see now that the general effect of an exponential distribution of sampling lifetimes is to 
decrease the peak width, and thus increase the coherency, by a factor of a/12 ks 3.5. 



In addition to the boxcar and Gaussian windows, another physically reasonable model 
for the hot spot evolution is that of a sharp rise followed by an exponential decay, perhaps 
caused by magnetic reconnection in the disk. In this case, the light curve would behave like 
a damped harmonic oscillator, for which the power spectrum is also given by a Lorentzian. 
Interestingly, the resulting QPO peak width is exactly the same, whether we use a collection 
of boxcar functions with an exponential lifetime distribution, or if we use a set of exponential 
sampling functions, each with the same decay time. In the discussion below, we will assume 
a boxcar sampling function and an exponential lifetime distribution, with its corresponding 
Lorentzian power spectrum. This also facilitates a direct comparison with observations, 
where the QPO data is often fit by a collection of Lorentzian peaks (Nowak 2000; Belloni et 
al. 2002). 

Due to the linear properties of the Fourier transform, the above analysis, while derived 
assuming a purely sinusoidal signal with a single frequency v , can be applied equally well to 
any periodic light curve with an arbitrary shape. If each coherent section of the light curve 
is written as 

f(t) = ^Ajsinpnvjt + </>]), (19) 

3 

then the total power spectrum (integrating over a distribution of coherent segments with 
random phase) is simply the sum of the Lorentz-broadened peaks: 

poo 

/» = / \G(u,T)\ 2 dN(T) 
Jo 

= 2/V, p ^|yj i + 4r2r ^_^ )2 . (20) 

Note that every peak in the power spectrum l 2 {y) has the same characteristic width Av = 
1/(2tt71). 
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The sum in equation (19) can be generalized to a Fourier integral so that (20) becomes 
the convolution of the segment power spectrum F 2 (u) with a normalized Lorentzian £(0, Au) 
centered on u = with width Au: 



Now we can apply our results to the light curves as calculated by the original ray- 
tracing code for a single geodesic hot spot. First, the X-ray signal over one complete period is 
calculated to give the Fourier components Aj in (19). For geodesic orbits, the power spectrum 
F 2 {u) is concentrated at integer combinations of the black hole coordinate frequencies u^, 
u r , and uq. In Schnittman & Bertschinger (2004a) we showed how the Fourier amplitudes 
depend on orbital inclination and eccentricity as well as hot spot shape and overbrightness. 
Given these frequencies Uj, amplitudes Aj, and a characteristic hot spot lifetime TJ, the 
integrated power spectrum follows directly from equation (20). 

Using the same ray-tracing code, we can also directly simulate the light curve and 
corresponding power spectrum produced by many hot spots orbiting with random phases, 
continually formed and destroyed over each time step with probability dt/Ti (Schnittman 
& Bertschinger 2004b). The power spectrum of such a simulation is shown in Figure 1 
(crosses), along with the analytic model (solid curve). We should stress that this curve 
is not a fit to the simulated data, but an independent result calculated using the model 
described above. For this particular example, the black hole has mass M = 1OM and spin 
a/M = 0.5 with a disk inclination of % — 70°. Each hot spot is on an orbit with u^ = 285 
Hz, u r = 95 Hz, and a moderate eccentricity of e = 0.1. This particular orbit was chosen 
because of the 3:1 ratio in coordinate frequencies, giving the strongest power at the modes 
z/0 — u r : z/0 : u^ + u r = 2 : 3 : 4. Similar frequencies appear to be the dominant peaks in the 
type A QPOs observed in XTE J1550-564 (Remillard et al. 2002). 

The defining characteristic of QPO peaks broadened by the summation of random phases 
is the uniform width of the individual peaks. For a power spectrum with multiple harmonics 
and beat modes, each peak is broadened by exactly the same amount, determined by the 
average lifetime of the individual hot spots. Thus if we can measure the widths of multiple 
QPO peaks in the data, the hot spot lifetime can be determined redundantly with a high 
level of confidence. 




(21) 
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3. DISTRIBUTION OF COORDINATE FREQUENCIES 

In the previous section, we assumed a single radius for the hot spot orbits. This assures 
identical geodesic coordinate frequencies for different hot spots with different phases and 
lifetimes. However, this assumption betrays one of the major weaknesses of the geodesic hot 
spot model: there still does not exist a strong physical argument for why these hot spots 
should form at one special radius or why that radius should have coordinate frequencies with 
integer commensurabilities. For now, we will be forced to leave that question unanswered, 
but we can make progress by drawing on intuition gained from other fields of physics. If 
there does exist some physical resonance in the system that favors these orbits, causing 
excess matter to "pile up" at certain radii (Abramowicz & Kluzniak 2001, 2003), then just 
like any other resonance, there should be some finite width in phase space over which the 
resonant behavior is important. The integer commensurability of the QPO peaks suggests 
that closed orbits may be playing an important role in the hot spot formation. If this is so, 
then some hot spots should also form along orbits that almost close, i.e. geodesies with nearly 
commensurate coordinate frequencies. These orbits will have guiding center radii similar to 
the critical radius r for which the geodesies form closed curves. 

Motivated by other processes in nature such as damped harmonic oscillators and atomic 
transitions, we model the resonance strength as a function of radius with a Lorentzian of 
characteristic width Ar. Then the probability of a hot spot forming at a given radius is 
proportional to the strength of the resonance there, giving a distribution of orbits according 
to 

P{r)dr= ^/frAr) 

For a relatively small resonance width Ar, we can linearize the coordinate frequencies Vj(r) 
around r = r with a simple Taylor expansion: 



u^r) « v jQ + (r - r ) 



(23) 



ro 



in which case the probability distribution in frequency space is also a Lorentzian: 



pw*,, = -w^ . ( 24 ) 

V i- V JQ 



Here Vj = v^, ug, v r are the azimuthal, vertical, and radial coordinate frequencies and Vj = 
Vj(ro) are those frequencies at the resonance center. 

For nearly circular orbits, the coordinate frequencies (using geometrized units with 
G — c — M — 1) are given by Merloni, Vietri, & Stella (1999) [following earlier work by 
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Bardeen, Press, & Teukolsky (1972); Perez et al. (1997)]: 



±1 



(25a) 



2vr(r 3 / 2 ± a) ' 



3a 



21 1/2 



v e = v 4> 1 =f 



(25b) 



and 



6 , 8a 




21 1/2 



v r = 1 



(25c) 



where the upper sign is taken for prograde orbits and the lower sign is taken for retrograde 
orbits (the results below assume prograde orbits, but the analysis for retrograde orbits is 
essentially the same). These frequencies are plotted in Figure 2 as a function of r for a 
representative black hole with mass 1OM and spin a/M = 0.5. The radial frequency 
approaches zero at the ISCO, where geodesies can orbit the black hole many times with 
steadily decreasing r. In the limit of zero spin and large r, the coordinate frequencies reduce 
to the degenerate Keplerian case with — v$ — v r — 1/ (27rr 3 / 2 ). 

Generally, the power spectrum of the periodic light curve from a single hot spot orbiting 
at r will be made up of delta functions located at the harmonics of the fundamental and 
the beat modes with v T and ug. Considering for the moment only planar orbits, the power 
will be concentrated at the frequencies v = nv^±v r , where n is some positive integer. In fact, 
there will be additional peaks at v — nz/^,±2z/ r and even higher beat-harmonic combinations, 
but for coordinate frequencies with = 3u r , these higher modes are degenerate, e. g. 
Vj, + 2u r = 2^ — v r . A careful treatment can distinguish between these degenerate modes, 
but in practice we find the power in the radial double- and triple-beats to be insignificant 
compared to the single-beat modes at nz/^ ± v r , so we limit our analysis to these frequencies. 

From equations (23) and (24), we see that a QPO peak centered around v = nu^ ± v r 
will be a Lorentzian of width 



Unlike in the previous section where the random phases gave a single width for every QPO 
peak, now each peak in the power spectrum will be broadened by a different but predictable 
amount. Note in particular how the peaks at the higher harmonics with n > 1 will be sig- 
nificantly broader (and thus lower in amplitude) than the fundamental. Another important 
feature evident from Figure 2 and equation (26) is that, due to the opposite-signed slopes of 
v r (r) and i^(r) around tq, the beat mode at v$ + v r remains very narrow, while the peak at 




(26) 
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— v r is quite broad. These features should play a key role in using the power spectrum as 
an observable in understanding the behavior of geodesic hot spots. 

As in Section 2, the first step in producing a simulated power spectrum is to calculate 
the Fourier amplitude in each mode with the full three-dimensional ray-tracing calculation 
of emission from a single periodic hot spot at r . Again, the linear properties of the problem 
allow us simply to sum a series of Lorentzians, each with a different amplitude, width, and 
location (peak frequency), to get the total power spectrum. The amplitudes Aj are given by 
the ray-tracing calculations, the locations isj from the geodesic coordinate frequencies and 
their harmonics, and the widths Auj from equation (26). 

Since the QPO peak broadening is most likely caused by a combination of factors in- 
cluding the hot spots' finite lifetimes as well as their finite radial distribution, the simulated 
power spectrum should incorporate both features in a single model. Now the computational 
convenience of Lorentzian peak profiles is clearly evident, since the net broadening is given 
by the convolution of both effects and the convolution of two Lorentzians is a Lorentzian: 

[C(u 1 ,Au 1 )*C(u 2 ,Au 2 )](u)=C(u tab ,Au tot )(u), (27) 

where the peak centers and widths simply add: u tot — v\ + v 2 and Au tot = Aui + Av 2 . In the 
case where one or both of the Lorentzians is not normalized, the amplitude of the convolved 
function is given as a function of the individual peak amplitudes and widths: 

A x A 2 Av x Av 2 

Aot = n Av l + Av 2 ' (28) 

where A\ and A 2 are the peak amplitudes of the respective Lorentzians [Aj = l/(jrAi/j) 
corresponds to a normalized function.] 

Figure 3 shows the simulated power spectrum for a collection of hot spots orbiting near 
the commensurate radius tq = 4.89M with a distribution width of Ar = 0.05M. All other 
black hole and orbital parameters are identical to those in Figure 1. Both the random phase 
broadening described in Section 2 and the effects of a finite resonance width are included in 
the model. Again, we should stress that the solid line is not a fit to the data, but rather an 
analytic model constructed from the sum of Lorentzian profiles as described above. In this 
example, the hot spots have a typical lifetime of 30 orbits, so the random phase broadening 
contributes only Av ~ 1.5 Hz for each peak. This allows us to focus on the effect that a 
finite resonance width has on the behavior of the QPO peaks at the coordinate frequencies 
and their various beat harmonics. For a resonance width of Ar = 0.05M, the peak widths 
due only to coordinate frequency broadening are shown in Table 1. 

The narrow peak at v$-\-v r = 380 Hz and the neighboring broad peak at 2z/^ — v r = 475 
Hz are clearly visible in the simulated data of Figure 3. Precise measurements of each 



- 12 - 



peak's width may not come until a next generation X-ray timing mission, but the qualitative 
behavior shown here should be detectable with the current observational capabilities of 
RXTE. Combining equations (11) and (26) gives a system of equations that solve for the 
hot spot lifetime TJ and the resonance width Ar as a function of the QPO peak widths 
Avj. If we could accurately measure the widths of only two peaks, both T\ and Ar could 
be determined with reasonable significance. More peaks would give tighter constraints and 
thus serve to either support or challenge the assumptions of the hot spot model. 



4. PHOTON STATISTICS AND CORONAL SCATTERING 

Another feature that is important in modeling QPO power spectra is based on the 
limitations of photon counting statistics by the RXTE instruments. We can easily model 
this effect using the same ray-tracing program employed above, now selecting only a small 
random fraction of the rays to make up the light curve at each time step. The resulting light 
curve is the "standard" light curve, effectively multiplied by a random function in time with 
amplitude described by a Poisson distribution, analogous to the sampling functions used in 
Section 2. The corresponding Fourier transform is the collection of original delta functions 
at the coordinate frequencies, convolved with the transform of a Poisson-valued sampling 
function. 

For a Poisson distributed function w(t) sampled over N intervals of At with an average 
value of fi, the corresponding Fourier transform W(v) is given by 

where for non-zero values of is, W(v) is a complex random variable with uniformly distributed 
phase and with amplitude given by the Rayleigh distribution [e. g., Groth (1975)] 

P(W)dW = —We- w2N/ »dW. (30) 
fi 

The corresponding power W 2 {y ^ 0) appears as a background white noise with mean 

(W 2 ) = £ (31) 

and variance 

AW 2 }=4^- 2 (32) 

constant over all frequencies. Unlike some of the other instrumental contributions to the 
background power, the Poisson noise is uncorrelated at different frequencies. 
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As mentioned in the previous sections, the transform of a perfectly periodic time series 
f(t) will be zero everywhere except for a finite set of frequencies Uj where F(uj) = A y To 
compute the simulated power spectrum of f(t), multiplied by a Poisson sampling function 
w(t), we must first convolve the transform functions F(y) and W(v) in frequency space 
and then take the squared amplitude, giving the net power spectrum G 2 (u), which will be 
composed of the same delta functions that define F 2 (u), scaled by a factor of ji 2 . In addition 
to these peaks, there will be a flat background noise function G 2 {y ^ Vj) with the same type 
of probability distribution as W 2 . The net power spectrum is then given by 

g2 m - { p { §) m 

where the noise spectrum G 2 {y ^ Uj) is again a random variable with distribution function 

NG 2 



P{G 2 )dG 2 oc exp ( - N ^ ) (34) 



and average power 

(G 2 \ = 

N 



(G 2 > = K4 ( 35 ) 



3 



One significant conclusion from this analysis is that the Poisson sampling noise will not 
contribute to the broadening of the QPO peaks. Since the Poisson noise is uncorrelated, 
a random sampling window function combined with any light curve [not just the strictly 
periodic f(t) used here] will maintain the Fourier properties of the original function F{y) 
while adding a constant level of background noise characterized by P(W). Thus we can 
confidently apply our analysis from Sections 2 and 3 without worrying about additional 
broadening from random emission and detection processes. Furthermore, these results give 
us insights into an equally important part of the QPO power spectrum: the background 
noise, which is critical when attempting to fit Lorentz functions to the data and determining 
the amplitude of each significant peak. 

From equations (33) and (35), we see that signal G 2 {y = Vj) to noise G 2 {y ^ Vj) scales 
as /i, thus giving a stronger detection for higher sampling rates, just as expected. In the 
context of the hot spot model, we interpret \x to be the average fraction of simulated light rays 
traced from the disk to the detector, so typically /i< 1. The actual value for fi will depend 
on the spatial and time resolution of the ray-tracing calculation as well as the luminosity 
and distance to the source in question. For a source like XTE Jl 550-564 radiating with an 
intensity of 1 Crab unit, the RXTE photon count rate in the range 2-60 keV is ~ 13, 000 
counts s _1 (Swank 1998). With a typical QPO amplitude of 1% rms, this corresponds to 
roughly 10~ 2 — 10° photons from the hot spot alone per 50 /is interval. 
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If the background noise were due purely to Poisson statistics, sufficient binning in the 
frequency domain would effectively remove the variance of G 2 and allow us to subtract the 
mean background, leaving a pure signal G 2 {vj). However, for the proportional counter array 
instrument on RXTE, there is an additional noise contribution from the effect of detector 
"deadtime" that must be considered during the data analysis process (van der Klis 1989; 
Mitsuda & Dotani 1989). This additional correction is particularly important for high count 
rates. For these high intensity sources, there are a number of other detector contributions 
that must also be included. In practice, they are generally combined and modeled as a 
single noise function introduced in the data analysis pipeline, such as a broken power law in 
frequency. 

Another simplified model we have included is that of scattering photons from the hot 
spot off of a low-density corona of hot electrons around the black hole and accretion disk. 
This is known to be an important process for just about every observed state of the black 
hole system (McClintock & Remillard 2004). Unfortunately, it is also an extremely difficult 
process to model accurately. Fortunately, for the problem of calculating light curves and 
power spectra, a detailed description of the scattering processes is probably not necessary. 
The most important qualitative feature of the coronal scattering is a smearing of the hot 
spot image: a relativistic emitter surrounded by a cloud of scattering electrons will appear 
blurred, just like a lighthouse shining its beam through dense fog. The effect is even more 
pronounced in the black hole case, where the hot spot orbital period is of the same order as 
the light-crossing time of a small corona, thus spreading out the X-ray signal in time as well 
as space. 

Since the scattered photons are often boosted to higher energies, a coherent phase lag 
in the light curves from different energy channels could be used to estimate the overall scale 
length of the corona. Vaughan et al. (1997) have observed this effect in neutron star QPOs 
and infer a scattering length of A ~ 5 — 15M for an optical depth of r ~ 5 in the source 
4U 1608-52. Ford et al. (1999) perform a similar analysis for black holes, including the 
possibility for an inhomogeneous corona, and derive a much larger upper limit (A ~ 10 3 M) 
for the size of the corona. In either case, the qualitative effect will be the same: the damping 
of higher harmonic features in the power spectrum of the X-ray light curve. 

The simple model we introduce is based on adding a random time delay to each photon 
detected from the hot spot. The distribution of this time delay is computed as follows: we fix 
the optical depth to be unity for scattering though a medium of constant electron density, 
so each photon is assumed to scatter exactly once between the emitter and the observer, 
thus determining the length scale of the corona as a function of density. In this case, the 
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probability of scattering after a distance / is 

P{l)dl = (36) 
A 

where A is the photon mean free path in the corona. 

Next, due to the likely existence of an optically thick disk around the black hole equator, 
we assume that the photon scattering angle is less than it/ 2 (we define the scattering angle 
9' as the angle between the incoming and outgoing wave vectors, so a straight path would 
correspond to 6' = 0). In other words, only photons emitted in a hemisphere facing the 
observer can ultimately be scattered in the observer's direction. For a photon emitted at an 
angle 9 to the observer, scattering at a distance / from the source produces an additional 
photon path length of d — 1(1 — cos 9), assuming for simplicity a flat spacetime geometry. 
While the photons are emitted with an isotropic distribution in, the scattering distribution 
is not isotropic. Since the scattering geometry requires that 9 = 9', we only detect a subset 
of the photons emitted with an angular distribution in 9 that satisfies this relationship. In 
the limit of low-energy photons (hv <C m e c 2 ) and elastic scattering, the classical Thomson 
cross section a T is used: 

^ = ^ T sin^(l + cos 2 ^)- (37) 

Integrating this distribution over all forward-scattered photons (9' < 7r/2), we find the 
average additional path length to be (d) = 7//16. Since the time delay is the path length 
divided by the speed of light c, scattering once in the corona adds a time delay At to each 
photon with probability 

P(At)d(At) = ^M e -At/T scat ? (3g) 
J- scat 

where the average scattering time is given by T scat = 7A/16c. 

Applied to the ray-tracing model, this has the effect of smoothing out the light curve 
with a simple convolution in the time domain of the original signal fit) and the time delay 
probability distribution function P(At). The Fourier transform of the resulting light curve 
is the product of the two transforms F(u) and -P(^), where 

P(y) = ^ — • ( 39 ) 

When we square the product to get the power spectrum G 2 (u) = F 2 (z/)P 2 (z/), the scaling 
factor is yet again a Lorentzian: 

^ - i + (40) 
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where the scale of frequency damping is given by 

Az/ scat = l/(27rT scat ) (41) 

and Aj are the delta function amplitudes of F{v) as defined above. This analytic result 
is perhaps a case where the ends justify the means. Our model for electron scatting in 
the corona is extraordinarily simplified, ignoring the important factors of photon energy, 
non-isotropic emission, multiple scattering events in a non-homogeneous medium, and all 
relativistic effects. However, assuming that almost any analytic model would be equally 
(in)accurate, at least the treatment we have applied proves to be computationally very 
convenient. 

Equation (40) states that the resulting power spectrum of the scattered light curve is a 
set of delta functions, with the higher harmonics damped out by the effective blurring of the 
hot spot beam propagating through the coronal electrons. A simulated power spectrum is 
shown in Figure 4a for a scattering length of A = 10M, comparable to the size of the hot spot 
orbit. Figure 4b shows the effect of a larger, low-density corona with scale length A = 100M, 
corresponding to a longer convolution time and thus stronger harmonic damping. The white 
background noise (Poisson noise with // = 1) in both cases is due to the statistics of the 
random scattering of each photon from one time bin to another. The simulated spectra are 
plotted as dots (asterices at Uj to highlight the peaks) and the analytic model is a solid line. 

As in the model for Poisson sampling, we see that the coronal scattering should not 
contribute to the broadening of the QPO peaks. However, it will have a very significant effect 
on the overall harmonic structure of the power spectrum, particularly at higher frequencies. 
Schnittman & Bertschinger (2004a) show a similar result caused by the stretching of the 
geodesic blob into an arc along its path, also damping out the power at higher harmonics. In 
this context, it is now clear that the arc damping can be modeled analytically by interpreting 
the stretching of the blob in space as a convolution of the light curve in time. If the stretched 
hot spot has a Gaussian distribution in azimuth with length A<f>, the original X-ray light 
curve will be convolved with a Gaussian window of characteristic time At = A0/(27rz^). 
Equation (13) gives the corresponding scaling factor W(v) in the frequency domain (replacing 
T with At/2). The exponential damping of W(v) is stronger than the Lorentzian factor at 
higher frequencies, but both effects (coronal scattering and hot spot stretching) are probably 
important in explaining the lack of significant power in the harmonics above ~ 500 Hz in 
the RXTE observations. 
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5. FITTING QPO DATA FROM XTE J1550-564 

In this section we combine all the pieces of the model developed above and apply the 
results to the RXTE data from type A and type B QPOs observed in the low-mass X-ray 
binary XTE J1550-564. To compare directly with the data from Remillard et al. (2002), we 
need to change slightly our normalization of the power spectrum. Following Leahy (1983) and 
van der Klis (1997), we define the power spectrum Q{v) so that the total power integrated 
over frequency gives the mean square of the discrete light curve Ij = I(tj): 



I 



1 ^ flj-(l) 



u>0 - o j=Q 

where Ij is sampled over j — 0, N s — 1 with average value (I). In terms of the power 
spectra used in Sections 2 and 3, Q{v) is given by 

Q(0) = 2 

Q{u) = 2T f C^ : (43) 
V ; ; 7 2 (0) V ; 



which has units of [(rms/mean) 2 Hz 1 ]. 

As we described in the introduction, the hot spot model is constructed in a number 
of steps. These steps result in a first approximation for the black hole and hot spot model 
parameters, after which a x 2 minimization is performed to give the best values for each data 
set. 



• The black hole mass and the inclination of the disk are given by optical radial velocity 
measurements (Orosz et al. 2002). We take M = 1O.5M and % = 72° as fixed in this 
analysis. 

• The black hole spin is determined by matching the frequencies of the HFQPOs to 
the geodesic coordinate frequencies such that v$ = 3u r at the hot spot orbit, giving 
a/M ss 0.5 for v$ sa 276 Hz. The small uncertainties in the measured value of can 
thus be interpreted indirectly as constraints on the mass-spin relationship. 

• The orbital eccentricity and hot spot size and overbrightness are chosen to match the 
total amplitude of the observed fluctuations. We use a moderate eccentricity of e = 0.1, 
but find the peak amplitudes are not very sensitive to this parameter (Schnittman & 
Bertschinger 2004b). The question of overbrightness is still an area of much research, 
since the nature of the background disk is not well known during the "steep power law" 
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state that produces the HFQPOS (McClintock & Remillard 2004). In practice, we set 
the hot spot emissivity constant and then fit an additional steady-state background 
flux I B to the variable light curve. 

• The hot spot arc length and the coronal scattering time scale are chosen to fit the 
relative amplitudes of the different QPO peaks. 

• The hot spot lifetime and the width of the resonance around tq are chosen to fit the 
width of the QPO peaks. 

• The observed flux gives the detector count rate, which is added to the white noise from 
the coronal scattering statistics. In practice, for power spectra with sufficient binning 
statistics (corresponding to long observations in time), the Poisson noise reduces to a 
flat spectrum which is subtracted from the data before fitting the model parameters. 

• As a final step, we include an additional power law component oc u^ 1 to account for the 
contribution due to turbulence in the disk [e. g. Mandelbrot (1999)] not accounted for 
by the hot spot model. Instrumental effects such as the detector deadtime are combined 
with the turbulent noise to give a simple two-component background spectrum: 

QnoiseH = QPL^" 1 + Qflat- (44) 

After determining the fixed amplitudes Aj with the ray-tracing calculation, we mini- 
mize the x 2 fit 

over the following parameters: orbital frequency z/^, hot spot lifetime TJ, 
resonance width Ar, scattering length A, hot spot arc length A0, steady state flux Ib, and 
the background noise components Qpl and Qflat- The best fit parameters are shown in Table 
2, along with la (68%) confidence limits. These confidence limits are determined by setting 
Ax 2 < 7.04, corresponding to six "interesting" parameters of the hot spot model, holding 
the noise components constant (Avni 1976; Press et al. 1997). We find that Qpl and Q flat 
are almost identical for both data sets, supporting the presumption that they are indeed a 
background feature independent of the hot spot model. 

In Figure 5 we show the observed power spectra for type A and type B QPOs, as 
reported in Remillard et al. (2002), along with our best fit models. The type A QPOs are 
characterized by a strong, relatively narrow peak at pa 280 Hz, corresponding to in 
our model, with a minor peak of comparable width at v$ — u r Pd 187 Hz. Type B QPOs 
on the other hand, have a strong, broad peak around 180 Hz with a minor peak at 270 Hz. 
This implies a longer arc, damping out the higher frequency modes, and a shorter average 
lifetime, broadening the peaks. Both types of QPO suggest a very narrow resonance width 
Ar, yet the current data does not constrain this parameter very well. Thus we assume the 
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majority of the peak broadening is caused by the addition of multiple hot spots with random 
phases and a characteristic lifetime of ~ 3 orbits for the type A QPOs and about half that 
for type B. 

We performed a covariance analysis of the parameter space near the y 2 minimum to 
identify the best-constrained parameters and their relative (in) dependence. This analysis 
confirms what the confidence limits suggest: the best-constrained parameters are the orbital 
frequency u^, the hot spot lifetime TJ, the arc length A(f>, and the background flux Ib- For 
the type A QPOs, we find and 7] to be independent, while the arc length and background 
flux are strongly correlated, so that A<P/Ib is positive and roughly constant within our 
quoted confidence region. This is because, for shorter arcs with fixed emissivity, increasing 
the arc length will increase the amplitude of the light curve modulation, requiring a larger 
background flux to give the same QPO amplitude. For the type B QPOs on the other 
hand, a longer arc length does not significantly amplify the modulation, since in the limit 
A0 — > 360°, the light curve would remain constant, and thus the parameters A0 and Ib 
are relatively independent. For both type A and type B QPOs, we find that the resonance 
width and the coronal scattering length are independent, yet not very well constrained. 

The resulting amplitudes and widths of the major QPO peaks are shown in Table 3, 
along with la confidence limits. These amplitudes are given by the analytic model so that 
the total rms in the peak centered at Vj is 

rm Sj = V2^-, (45) 

where A' Q is the mean amplitude of the light curve (including the background Ib) and A'- 
are the original Fourier amplitudes Aj given by the ray-tracing code, appropriately scaled 
according to equations (33) and (40). This is more instructive than measuring the rms 
directly from Q(z/), which includes background power and instrumental effects uncorrelated 
to the actual QPO peaks. 

In Schnittman & Bertschinger (2004a), the hot spot light curve was added to a steady- 
state disk with emissivity that scales as r~ 2 , which provides an estimate of the size and 
overbrightness of the hot spots required to produce a given (rms/mean) amplitude in the 
light curve. Considering that most high frequency QPOs are observed with the greatest 
significance in the 6-30 keV energy band during the steep power-law spectral state (McClin- 
tock & Remillard 2004), it seems rather unlikely that the background flux is coming from 
a thermal, optically thick disk. Even if the flux is originally produced by such a disk, it 
clearly undergoes significant scattering in a hot corona to give the high temperature power 
law observed in the photon energy spectrum. 

In the context of the model presented here, we can only calculate the fraction of the 
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total flux that is coming from the hot spots, determined by fitting to the QPO data, without 
presuming an actual model for the background emission. For XTE J 1550-564, we find that 
the type A hot spot/arcs contribute 8.5% of the flux in the 6-30 keV band, while the type 
B arcs must contribute significantly more (38%) to give a comparable amplitude. This is 
due to the longer arc length described above: in the limit of an azimuthally symmetric ring, 
even infinite brightness would produce no variability. 

6. DISCUSSION AND CONCLUSIONS 

In the context of a geodesic hot spot model, we have developed a few simple analytic 
methods to interpret the amplitudes and widths of QPO peaks in accreting black holes. 
The model combines three-dimensional ray tracing calculations in full general relativity with 
analytic results of basic convolution theory, which are in turn confirmed by simulating the 
observed light curves of multiple hot spots. Given the Fourier amplitudes of a single hot 
spot, we have derived a simple formula for the complete QPO power spectrum made up of 
Lorentzian peaks of varying amplitudes and widths. This power spectrum can then be fit to 
observed QPO data and used to constrain parameters of the hot spot model, and possibly 
measure the black hole mass and spin. 

For XTE J1550-564, the locations of the HFQPO peaks are well constrained, in turn 
constraining the spin parameter a/M when combined with radial velocity measurements 
of the black hole mass. Based on the presumption that the 3:2 frequency ratio is indeed 
caused by closed orbits with coordinate frequencies in a 3:1 ratio, an observed mass of M — 
10.5 ± l.OM and orbital frequency = 276 ±5 Hz would predict a spin of a/M = 0.5 ±0.1 
(Orosz et al. 2002; Remillard et al. 2002). If reliable, this coordinate frequency method 
would give one of the best estimates yet for a black hole spin. 

The amplitudes of the QPO peaks can be used to infer the arc length of the sheared hot 
spot and the relative flux contributions from the hot spot and the background disk/corona. 
The longer arcs seen in type B QPOs are also consistent with the broader peaks: if the 
hot spots are continually formed and destroyed along special closed orbits, as the emission 
region gets stretched into a ring, it is more likely to be dissolved or disrupted, giving a shorter 
characteristic lifetime T t and thus broader peaks. 

Unfortunately, the quality of the QPO data is not sufficiently high to confirm or rule out 
the present hot spot model, leaving a number of questions unanswered. By fitting only two 
or three peaks, we are not able to tightly constrain all the model parameters, particularly 
the scattering length scale and the resonance width, both of which are most sensitive to 
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the higher frequency harmonics. Since the high temperature electrons in the corona tend 
to transfer energy into the scattered photons, measuring the energy spectra of the different 
QPO peaks would also prove extremely valuable in understanding the emission and scattering 
mechanisms. This has been done to some degree with the lower frequency region of the power 
spectra from black holes and neutron stars (Ford et al. 1999), and may even be observable 
above ~ 100 Hz with current RXTE capabilities. For this analysis to be most effective, a 
more accurate model for the electron scattering will certainly be necessary. 

Some of the power spectrum features discussed in this paper are unique to the geodesic 
hot spot model, while others could be applied to more general QPO models. Clearly the 
harmonic amplitudes Aj given by the ray-tracing calculation are dependent on the hot spot 
model, as is the broadening from a finite resonance width, yet both could be generalized 
and applied to virtually any perturbed disk model. Similarly, the random phase broadening 
and the damping of higher harmonics due to photon scattering will be important effects for 
any emission mechanism that produces periodic light curves from black holes. If the next 
generation X-ray timing mission could produce power spectra comparable to the phenome- 
nal cosmic microwave background (CMB) results of recent years, we believe that many of 
these issues could be resolved. Energy resolution and polarization would similarly provide 
extremely valuable information about the source of the QPOs. As with the CMB, each 
successive peak of the power spectrum would help to pin down another parameter until the 
model becomes predictive instead of descriptive, or is ruled out all together. 

In the immediate future, however, there is much more to be done with the RXTE data 
that already exists. Important additional insight might be gained from new analysis of the 
X-ray light curves in the time domain, recovering some of the phase information lost when 
the power spectrum is computed in frequency space. There is also an important message 
in the relationship between the photon energy spectra and the QPO power spectra as well 
as the connection between the low frequency and high frequency QPOs. Why should the 
HFQPOs appear in certain spectral states and not others? The answer to these questions 
may lie in new models of the accretion disk and specifically the radiation physics relating 
the thermal and power-law emission, as well as broad fluorescent lines like Fe Ka. The fact 
that the HFQPOs are seen most clearly in the 6-30 keV energy range suggests that standard 
models of thin, thermal accretion disks are not adequate for this problem. This emphasizes 
the essential role of radiation transport, particularly through the corona, in any physical 
model for black hole QPOs. 



We thank Ron Remillard for many helpful discussions and providing the QPO data for 
XTE J1550-564. This work was supported by NASA grant NAG5-13306. 
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Table 1: Widths of QPO peaks around coordinate frequency modes nv^ ± z/ r , due to a radial 
distribution of hot spots with Ar = 0.05M. For relatively narrow resonance regions, the 
QPO peak widths are linearly proportional to Ar. The basic black hole and hot spot model 
parameters are the same as in Figures 1 and 3. 

Mode Frequency (Hz) FWHM (Hz) 



v r 95 3.6 

- v r 190 12.2 

i/0 285 8.4 

v$ + v r 380 4.8 

2z/^ - v r 475 20.6 

2z/^ 570 16.8 

2z/^ + v r 665 13.2 



Table 2: Best-fit parameters of the hot spot model for type A and type B QPOs from XTE 
J1550-564. (la) confidences are shown in parentheses. 



Parameter 


Type A 


TypeB 


orbital frequency (Hz) 


280.1(2.4) 


270.5(12) 


lifetime T/ (ms) 


10(2.0) 


5(1.5) 


(orbits) 


2.8(0.55) 


1.4(0.4) 


resonance width Ar (M) 


0.02(0.05) 


0.025(0.12) 


scattering length A (M) 


5(10) 


10(20) 


arc length Ac/> (°) 


155(30) 


285(20) 


flux ratio T 7 !;° tspot 

JS+Jhotspot 


0.085(0.025) 


0.38(0.05) 
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Table 3: Amplitudes and widths of type A and type B QPO peaks from XTE J1550-564, as 
determined by the best fit parameters listed in Table 2 and equation (45). (la) confidences 
are shown in parentheses. 

A B 



Mode 


rms 


FWHM 


rms 


FWHM 




(%) 


(Hz) 


(%) 


(Hz) 


v r 


0.57(0.15) 


33.1(6.2) 


2.03(0.21) 


63.6(16.0) 


v<p - v r 


1.62(0.26) 


35.7(5.9) 


2.57(0.14) 


67.6(15.5) 


v<t> 


3.35(0.17) 


34.6(5.5) 


1.48(0.24) 


65.9(15.3) 


Z/0 + v r 


0.75(0.19) 


33.4(5.8) 


0.06(0.02) 


64.1(15.8) 
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Fig. 1. — Simulated power density spectrum (crosses) from a ray-tracing calculation of 
many hot spots on geodesic orbits with random phases and different lifetimes, along with an 
analytic model (solid line) of that power spectrum. The black hole has mass M = 1OM and 
spin a/M = 0.5, giving v r = 95 Hz and v$ = 285 Hz. The hot spot orbit has an eccentricity 
of 0.1 around a radius of r = 4.89M and an inclination of 70°. The peaks have Lorentzian 
profiles with Av « 11 Hz, corresponding to a characteristic hot spot lifetime of four orbits. 
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Fig. 2. — Geodesic coordinate frequencies as a function of radius for a black hole with mass 
M = 1OM and spin a/M = 0.5. The radius of the inner-most stable circular orbit risco 
is where v r — > 0. The commensurate radius r is where the ratio of azimuthal to radial 
coordinate frequencies is 3:1. 
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Fig. 3. — Power density spectrum (crosses) from a ray-tracing calculation of many hot spots 
on geodesic orbits with different radii r and thus different coordinate frequencies, along with 
an analytic model (line) of that power spectrum. The black hole has mass M = 1OM and 
spin a/M — 0.5, while the average hot spot orbit has an eccentricity of 0.1 around a radius of 
r = 4.89M, as in Figure 1. The peaks have Lorentzian profiles with Av given by equations 
(11) and (26) with T x = 100 ms and Ar = 0.05M. 
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Fig. 4. — Power density spectrum from a single hot spot light curve where the emitted 
photons are scattered exactly once each by a uniform corona of electrons. The simulated 
spectra are plotted as dots and asterices, while the analytic model is a solid line. In (a), the 
mean free path to scattering is A = 10M, while (b) represents a much larger, low density 
corona with A = 100M. 
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Fig. 5. — Comparison of hot spot model power spectrum (line) with data (crosses) from 
XTE J1550-564. (a) Type A QPO, dominated by a narrow peak at v$ « 280 Hz. (b) Type 
B QPO, dominated by a broad peak at v^ — v T ^ 180 Hz. The best fit model parameters for 
each data set are shown in Table 2 and the resulting QPO amplitudes and widths are shown 
in Table 3. 
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